Chiral Polaron Formation in Graphene 
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A theoretical investigation ol the possible existence of the chiral polaron formation in graphene 
is reported. We present an analytical method to calculate the ground-state of the electron-phonon 
system within the framework of the Lee-Low and Pines theory. On the basis of our model, the 
influence of electron-optical phonon interaction onto the graphene electronic spectrum is examined. 
In this paper, we only considered doubly degenerate optical phonon modes of E2 9 symmetry near 
the zone center F. We show analytically that the energy dispersions of both valance and conduction 
bands of the pristine graphene differ significantly than those obtained through the standard electron 
self energy calculations due to the electron-phonon interactions. Furthermore, we prove that the 
degenerate band structure of the graphene promote the chiral polaron formation. 
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Since the discovery of grapheneA, many studies have 
focused on electron-phonon interactions, because of their 
particular importance in understanding electronic and 
optical properties of graphene^—. Even though the 
interaction of electron with doubly degenerate optical 
phonon modes of Fi2 g symmetry near the zone cen- 
ter r does not open a gap^I, they contribute signifi- 
cantly to intravalley-intraband and intravalley-interband 
scattering™ Since these phonon modes consist of in- 
plane out-of-phase displacements of two sublattices A 
and B, they just shift the K (K ) point so as to electronic 
bands are still described by Dirac cones 4 . Moreover, it is 
also revealed that the gate-modulated low-temperature 
Raman spectroscopy, the graphene G band, which is the 
optical phonon at long wavelength, is markedly sensitive 
to the coupling with Dirac fermion excitations at small 
wave vectors (long wavelengths) 5 . 

So far only pcrturbational methods for intraband tran- 
sitions have been applied to the electron-phonon interac- 
tions, the combined effect of intra- and inter-band scat- 
tering in a single valley has not been considered. Due 
to the fact that graphene is a semimetal, actually gap- 
less or zero-gap semiconductor, the Dirac cones can not 
be treated independently for electron-E2 9 phonon inter- 
actions. Single band approximations take into account 
the intraband transitions only and they ignore the lack 
of gap, thus omit the electron (hole) transitions from it 
(it* ) valance (conduction) to ir* (tt ) conduction (valance) 
bands. Whereas, due to the lack of gap, F,2 g phonon can 
easily excite electron-hole pairs leading to the chirality 
dependent modifications in carrier dynamics. 

In this letter, we used the continuum Frohlich type 
model to treat the interaction of electron with long 
wavelength E 2ff phonon. We introduce a diagonaliza- 
tion procedure based on Lee-Low-Pines (LLP) like trans- 
formations to investigate the properties of both valance 
and conduction band polarons in graphene. Within the 
framework of low-energy continuum model of graphene, 
the Hamiltonian of an electron interacting with optical 
T-phonon around the K point in the Brillouin zone can 
be written as 



(1) 



where Hq — vper ■ p is the unperturbed bare Hamilto- 
nian, whose spectrum describes cone like behavior with 
eigenvalues ek\ — Ae^ wherein €k = vpk and A is the chi- 
rality index, together with the corresponding eigenkets 

|feA) = cxp(ifc • r) ( Xe i0 ^ ) f /y/2L. 

In Eq.(p}, H e -p is the electron-optical phonon interac- 
tion HamiltonianX i 13 ! 24 , and is given by 



E 



H.c. 



(2) 



where b qil (b qfl ) is the optical phonon creation (annihi- 
lation) operator with longitudinal and transverse opti- 
cal phonon branch index /i = 1 (LO) and 2 (TO), re- 
spectively. Their dispersion have the form tu^ (q) — 
ur (0) (qa) with dimensionless part (qa) — 

1/2 

[1 — r M (qa)] ' , and with wr (0) = 0.196 eV, where n 
and r2 are given by 0.1138 and 0.01979, respectively^. 
The momentum dependent matrix element of the inter- 
action Eq.© is (q) = (q) M M (q) , and defined 



M„ (q) 



M B Ae" #(q) 
M AB e^( g ) 



(3) 



with M AB = +l(i) and M BA = -1(«) for LO(TO) 
phonons--, <f> (q) = tan -1 (q y /q x ) is the azimuthal 

angle of the phonon wave vector q, and M„ (q) — 
M/^/NRf, (qa) with M = 3a q Jo/V^- Here, q = 
(dJo/da) /J is predicted2L2S around 2 I" 1 or 2.5 A' 1 . 

To solve Eq.([T]), we propose a diagonalization proce- 
dure based on the LLP method, which includes two suc- 
cessive unitary transformations. Pristine graphene has 
a electron-hole degeneracy point at k = 0. Therefore, 
to be compatible with this gapless band structure of 
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the graphene, we make an ansatz for the chiral polaron 
ground-state vector 



pol 



= U 1 U 2 \0) 



ph 



l±) 



such that H\$) 



I pol - ^± 

for the phonon vacuum 



and 



In Eq.©, |0) , stands 



±>=E; 



ph 

\kX) 



sponds to electronic state vector defined through the ap- 
propriate fractional amplitudes, a±, due to the fact that 
polaronic wavefunction must be the linear combination 
of |fe+) and |fe— ), respectively. While the first unitary 
transformation 



U\ = exp 



-ir 



eliminates the electron coordinates, since the trans- 
formed operators are given by the relations, b qfJi = 
b qiJ , exp [—iq ■ r] and p — p — J2 q ^>\j>q^ second uni- 



( A \ 

W tary transformation 



U 2 



exp 




K («) 



|feA) - H.c. 



is the displaced oscillator transformation with amplitude 


Mo M (9) = (q)/huj f j i (q). It just shifts the phonon 
coordinates. As a result, the transformed Hamilto- 
nian can be written as % = U2 1 U^ 1 T-LU\U2- Therefore, 
p h(0\rl \0) h leads to the following equation: 



ph 



(0\H\0) ph =v F cr- (p-hJ2q\Mon(q) 
\ q 

+ ^7kJ M (qr)|M 0M (q)| 



feA 



Mt(q)|fcA) -J2 M fl (q)\M 0lM (q) 



kX 



H.c. 
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Finally, by using the ansatz given by Eq.(j4]), one can easily construct characteristic equation of the matrix in 

the form 

I 



E± - hv F k + Sw F S ++ (^, k) + T}% + (ji, k) ftu F 53 + _(^, k) + 2£ t [_( ) u, k) 

Uv F 'S-+(n, k) + 2£°_ + (>, k) E± + hv F k + %v F Ji—(n, k) + E° _(^, k) 



(5) 



with elements 



£++0,A:) 
S+-(/i,fc) J 2 



JEI M °m(9)| 



x [l±s„ cos 2(6-0)] 



<?,i cos ( 
iq sin ( 



x hupfa) [l + s^XX cos 2{0-4>) , (6) 



and S_ + = S! 



After convert- 



where X = — S++ 

ing the sums in Eq.(j6|) into integrals over q, it is easy to 
see that the terms with (ji, k) do not give contribu- 
tion to the eigenvalue calculation, since they all vanish 
after the q integration. But the rest, i.e., terms with 
£^ A , (jU, k) contribute. Thus, from the secular equation, 
i.e., from Eq.© the eigenvalues E± can then be solved 
analytically as a function of k in closed form, 



E± = ± { {hv F kf + 4 [S° ±T (n, k)] 2 } 



1/2 



(7) 

In Eq.©, while the last term is due to the intraband 
transitions, second one in the parenthesis is contribu- 
tion due to the coupling between valance and conduc- 
tion bands, i.e. it corresponds to interband transitions. 
Thus, our chiral polaron dispersion consists of partial 
mixture of contributions from valance and conduction 
bands. While taking the integrals in Eq.JTJ), since they 
diverge at upper limit of the integrations, we must in- 
troduce an upper cut-off frequency in the integrations to 
make them to be finite. It is nothing but just 2k, due to 
q £ [0,2k]. As a result, one obtains k— dependent contri- 
butions as 



S AA'(^) = ^oa(Q)ln 



1 



vr \2 



(8) 



where we have defined a (0) = |M| /4Johu>r(0) and new 
dimensionless wave vector, k = ka. a (0) is of order 
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FIG. 1. (Color online) The polaron ground-state energy ob- 
tained from Eq.jTJ as a function of ko = ka. Inset: shows the 
electron dispersion normalized to Jo as a function of ko = ka. 



of 0.054 (0.089) depending on the choice of q = 2 (2.5) 
A _1 .We analyze Eq.0 in FIG.[TJ which shows both po- 
laron self energy and Fermi velocity renormalization. As 
it can be easily seen from the figure that polaron self en- 
ergy is strongly enhanced with increase in k Q . This can 
also be justified from its inset that shows how electron- 
phonon interactions renormalize the Dirac velocity. To 
see this clearly we also present in FIG. [2] electron-hole 
band dispersions as a function of ko — ka both in the 
absence and presence of electron (hole)- phonon interac- 
tions. 

It is convenient to rewrite Eq.© for small k's, which 
corresponds to neglect phonon dispersions, i.e., it results 
with polaron dispersion independent on r M . To do this, 
we first expand the logarithm in Eq.©, in power series 
of k , and then replace the resultant back into Eq.([S])so 
that Eq.0 reduces to the simple form 



E± = ±k a ( 1 + 4 



2V3 



a (0) k c 



1/2 



2 ^ /nil 2 
7T 



(9) 

where E± = aE±/frvp is the dimensionless energy. The 
standard electron-self energy calculations due to the in- 



teractions of electron with degenerate optical phonon 
modes with E 2s symmetry in graphene predict a lower- 
ing of both conduction and valance band energies in the 
same direction to preserve the symmetry of the Dirac 
cones. Whereas, in our case, as k Q increases, negativ- 
ity due to intraband interactions is compensated by the 
second term in the parenthesis in Eq.©, i.e., by the in- 
terband interactions. Although it is compensated, their 
combined effects strongly modify the Dirac cones (FIG. 
2 ). Since Eq.© can easily be rearranged into the form 
E\ = \fivp (ko) ko , we can thus define the renormalized 
Fermi velocity as 



vf (k ) = v F 



^Jl + Aa 2 (ko)Ta(k ) 



wherein, in analogy with quantum chromodynamics— , 
we also define a new coupling constant, running coupling 
constant, as a function of k, i.e, energy, 



E/./„ 




FIG. 2. Electron-hole band dispersions as a function of 
ko = ka both in the absence of (inner cone) and in the pres- 
ence of electron (hole)- phonon interaction (outer cone). Only 
LO phonon contribution is taken into account. The thick dark 
circle on the E/Jq — plane is the projection of both per- 
turbed (lower contour in the upper cone) and unperturbed 
(upper contour in the upper cone) the contour at E/ Jo = 1. 
It is drawn to guide to eye how the dressed polaronic band 
dispersions differ from the undressed ones. 



4 



a (fc ) = — a (0) fc - 

In conclusion, we predict a new type of polaronic 
formation in pristine graphene. We demonstrate that 
our results are differ from those obtained through stan- 
dard electron self energy calculations due to electron-E2 9 
phonon interactions in nondegenerate band case. We 
show that chiral polaron band dispersions consist of k 
dependent terms besides the free undressed one. In addi- 
tion to free undressed ones, both intraband and interband 
interactions coexist in a one polaron dispersion. More- 
over, a considerable renormalization of Fermi velocity, 
i.e., slope of Dirac cones in graphene is observed due to 



the electron-optical phonon interactions in graphene. It 
is also found that, the effect of LO phonon-electron in- 
teraction is stronger than that of TO. As for the validity 
of our approximation it is valid for fco = fca < 1, however 
it can be extended to higher fc values by just including 
the effect of trigonal warping in the total Hamiltonian. It 
can also be extended to the case electron-K phonons in- 
teraction where a gap occurs. Furthermore, it can also be 
generalized to the calculation of phonon induced electron- 
hole or electron-electron interactions, i.e., to exciton or 
bipolaron binding energies in graphene. 
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